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Abstract 

We consider the bosonic 5/(2) sector of the maximally supersymmetric ^ = 4 SYM 
model and show that anomalous dimension of the twist-3 single-trace composite op- 
erators built of scalar fields, recently calculated up to the four-loop order, can be gen- 
erated by a compact reciprocity respecting evolution kernel. 

1 Introduction 

QCD shares the vector boson — gluon — sector with supersymmetric Yang-Mills 
models (SYM). This suggests to explore supersymmetric partners of QCD in order 
to shed light on the subtle structure of the perturbative quark-gluon dynamics near 
the light cone. The latter manifests itself in parton distributions (structure functions 
of deep inelastic lepton-hadron scattering, DIS) and parton— >^hadron fragmentation 
functions (internal structure of jets in, e.g., e+e" annihilation into hadrons) which 
evolve with the hardness scale of the process, Q^. 

QCD is not an integrable quantum field theory. In spite of this, in certain sectors 
of the chromodynamics the integrability does emerge [1]. This happens, markedly, 
in the problem of high energy Regge behaviour of scattering amplitudes in the large 
Nc approximation (planar 't Hooft limit), in the spin | baryon wave function, for the 
scale dependence of specific (maximal helicity) quasi-partonic operators (for review 
see [2]). What all these problems have in common, is the irrelevance of quark degrees 
of freedom and the dominance of the classical part of gluon dynamics, in the sense of 
the Low-Burnett-Kroll theorem [3]. 

A SYM dynamics exhibits the deeper integrability, the higher the symmetry. The 
maximally supersymmetric jV = 4 YM theory occupies an exceptional position in this 
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respect. A string of recent impressive theoretical developments [2, 4, 5] hints at an 
intriguing possibility that this QFT, super-conformally invariant at the quantum level, 
may admit an all loop solution for anomalous dimensions of its composite operators. 
QCD would benefit a lot from such a solution, since this would provide a one-line- 
all-loops representation for the dominant part of the perturbative gluon dynamics, the 
part that could be dubbed LBK-classical. 

For the non-compact s[(2) sector containing twist operators with arbitrarily high 
spin the three-loop Bethe Ansatz equations (BAE) have been conjectured in [6]. The 
perturbative expansion of the Staudacher BAE perfectly matches the three-loop anoma- 
lous dimensions of twist-2 Wilson operators that determine the Bjorken scaling viola- 
tion pattern. At the same time, it apparently fails at the fourth loop where the wrapping 
problem is manifest [7]. 

Since in jV/" = 4 SYM all twist-2 operators belong to one super-multiplet, diago- 
nalization of the matrices 7(3 x 3) and 7(2 x 2) describing evolution of unpolarized 
(g, A, (p) and polarized parton distributions (g, A) results in the relation [8] 

j+iN + 2) = %{N + 1) = 7o(iV) = 7_(i\r - 1) = 7_(i\r - 2) = 7uni(iV). (1.1) 

It expresses all five diagonal elements in terms of the unique function — the "univer- 
sal anomalous dimension" 7urii with shifted arguments. Inspired by the structure of 
the answer in the first two loops [9-12], Kotikov, Lipatov, Onishchenko and Velizhanin 
(KLOV) have proposed in [13] the maximum transcendentality principle according to 

(n) 

which 7^- at n loops is a linear combination Euler-Zagier harmonic sums of transce- 

(31 

dentality r = 2n — 1. This allowed them to predict 7^^^^ by simply picking up from 
the three-loop non-singlet QCD anomalous dimension [14] the "most transcendental" 
terms (r = 5), having set Cp = Ca = N'c- 

(n) 

The complexity of higher loop expressions for 7^- , n > 2, is significantly reduced 
if one trades the anomalous dimension for a new object — the reciprocity respecting 
(RR) evolution kernel, [15]. In the case of 7^ = 4 SYM {P{a) = 0) this trade-off takes an 
exceptionally simple form. 



j,{N)=P{N-laj^{N)), a=\^\' , (1.2) 




where 7±(iV) are the anomalous dimensions responsible for time- and space-like par- 
ton evolution, respectively (in the S case which admits the operator product expan- 
sion, 7_ describes the scale dependence of Wilson operators of the leading twist 2). 
Here we employ the definition of the anomalous dimension which is more familiar to 
the integrable communit}0; within this convention the scaling violation rate is given 
by the expression 

^l^\nDM{Q') = -2^{N), (1.3) 



^The r.h.s. of (|1.31 is traditionally called 7 in the QCD parton evolution context. 
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with Dj^{Q^) the Mellin moment of the parton distribution measured at the large 
momentum transfer (hardness) scale Q^. 

In the context of the QCD parton picture, the notion of the RR evolution kernel V 
emerges as a result of the reformulation of space-like (DIS, "S") and time-like parton 
multiplication processes (e+e", "T") in terms of a unified evolution equation [16, 17]. 
This equation is constructed on a basis of the parton fluctuation time ordering. In the 
space of Bjorken (S)/Feynman (T) variable x conjugate to the Lorentz spin N, 

V{N) = x^V{x) = M [P{x)\ , (1.4) 

its integral kernel V{x), identical for the two channels, satisfies the Gribov-Lipatov 
reciprocity (GLR) [18] in all orders: 

F{x) = -xF{x-^). (1.5) 

Existence of the common evolution kernel V and its internal symmetry is reflected in 
the structure of the anomalous dimensions 7^ that one finds solving (|1.2|) . 

Basso and Korchemsky in [19] approached this problem from the point of view 
of the large N expansion. They have generalised ( 11.21 ) to anomalous dimensions of 
quasi-partonic operators of arbitrary twist L, and traced its origin to the conformal 
symmetry. 

Having analysed all anomalous dimensions of twist-2 operators, known to the 
two- and/ or three-loop order in QCD and SYM theories, as well as in the scalar X(p^ 
QFT (known at four loops), the authors demonstrated that the asymptotic series for 
the corresponding kernels V{N) run in integer negative powers of N{N + 1) — the 
quadratic Casimir of the collinear 5L(2; M) group. As a result, even, iV^^", and odd 
terms, N~^'^~^, of the large N expansion of anomalous dimensions 7o- turn out to be 
related. Basso and Korchemsky named this feature "parity preserving asymptotic 
series". 

In = 4 SYM the property of parity preserving series directly follows from the 
representation of the twist-2 kernel developed in three loops in [15]: 

PiN) = ASiiN)- (^^^^ + AiN)j +BiN). (1.6) 

Here aph is the physical coupling which determines the strength of LBK-classical 
radiation [16, 20] and coincides with the so-called cusp anomalous dimension [21] 
whose all-order weak-coupling expansion in 7\A = 4 SYM has been remarkably guessed 
P2]. 

The functions A = O (a^) and B = O (a^) have compact expressions in terms of 
RR combinations of complementary nested harmonic sums. Their asymptotic expan- 
sion at iV -> +CO is regular, that is, contains no In^ N enhancement factors at any level 
of the 1/N suppression. This makes the harmonic function 5i the only source of the 
log N behaviour in (|1.6|) , at least at the three-loop level. 
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This feature of the evolution kernel ( I1.6I ) is in a marked contrast with the anoma- 
lous dimension per se, whose large N expansion includes growing powers of log N: 

OO H fc 

j{N) = alnN + Y,j^Yl ^^rn In'" N. (1.7) 

fc=0 m=0 

Physically, the reduction of singularity of the large N expansion is due to the fact 
that towers of subleading logarithmic singularities in the anomalous dimension are 
actually inherited from the first loop — the LBK-classical ^'^^^ = V^^^ oc Si, and the RR 
evolution equation (|1.2|l generates them automaticall}|l [15, 17]. 

In this letter we consider twist-3 operators in the s[(2) sector of 7\A = 4 SYM and 
construct the evolution kernel V that generates the minimal anomalous dimension of 
single trace operators built of three scalar fields. This evolution kernel V satisfies the 
Gribov-Lipatov reciprocity, with substituted for x in eq. {1.5) , in accord with the 
fact that the harmonic functions entering the twist-3 anomalous dimension have JV/2 
for the argument [7,23]. We also demonstrate that the twist-3 kernel V admits the 
representation {1.6} , with A = O (a'*). 

This shows that the reciprocity respecting evolution, and thus the property of par- 
ity preserving series, manifest themselves also beyond the leading twist. 

2 Reciprocity Respecting evolution kernel 

2.1 Anomalous dimensions of twist-3 operators 

The sl(2) sector of planar jV = 4 SYM contains single trace states which are linear 
combinations of the basic operators 

Tr {(P^i Z) ■■■ {V^ Z)}, si + --- + SL = N, (2.8) 

where Z is one of the three complex scalar fields and D is a light-cone covariant 
derivative. The numbers {si} are non-negative integers and N is the total spin. The 
number L of Z fields is the twist of the operator, i.e. the classical dimension minus 
spin. The subsector of states with fixed spin and twist is perturbatively closed under 
renormalization mixing. 

The anomalous dimensions of states {2.8} are the eigenvalues ^L{N]g) of the di- 
latation operator — integrable Hamiltonian. These values were obtained by solving 
numerically the Bethe Ansatz equations (BAE), order by order in g^, and guessing the 
answer in terms of harmonic sums of transcedentality r = 2 n — 1, with n the number 
of loops. Since wrapping problems, delayed by supersymmetry, appear at L-l-2 loop 
order for twist-L operators [4, 24], the Bethe Ansatz equations for twist-3 are reliable 
up to four loops (including, at the fourth loop, the dressing factor). 

or the space-like case under consideration, cr = — 1 
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One finds [7,23] 

7?^ = 4 5i, (2.9a) 

7^'^ = -2(53 + 2 5i52) (2.9b) 
7^'^ = 5 55 + 6 52 53 - 8 53,i,i+4 54,i-4 52,3 + 5i(4 5| + 2 54 + 8 53,i), (2.9c) 
^ St + 7 Sifi + 15 52,5 — 5 53,4 — 29 54,3 — 21 55,2 — 5 56, 1 
-40 5i,i,5 - 32 5i,2,4 + 24 5i,3,3 + 32 5i,4,2 - 32 52,i,4 + 20 52,2,3 
+40 52,3,2 + 4 52,4,1 + 24 53,1,3 + 44 53,2,2 + 24 53,3,1 + 36 54,1,2 

+36 54,2,1 + 24 55,1,1 + 80 5i,i,i,4 - 16 5i,i,3,2 + 32 5i,i,4,i (2.9d) 



(4) 
73 



24 5i,2,2,2 


+ 16 5i,2,3,l 


- 24 5i,3,i,2 


-24 5i,3,2,i 


- 24 5i,4,i,i 


24 52,1,2,2 


+ 16 52,1,3,1 


— 24 52,2,1,2 


— 24 52,2,2,1 


- 24 52,3,1,1 


24 53,1,1,2 


- 24 53,1,2,1 


- 24 53,2,1,1 


-24 54,1,1,1 


- 64 5i, 1,1,3,1 


8^5i53. 











The last term in (|2.9d|) , with P = ^3, is the contribution from the dressing factor that 
appears in the BAE at the fourth loop [22]. 

The twist-3 anomalous dimension has two characteristic features: 

1. All harmonic functions 5a are evaluated at half the spin, 5a = 5a {N/2). On the 
integrability side, this does not look unwarranted, since only even N belong to 
the non-degenerate ground state of the magnet. 

2. No negative indices appear in (|2.9|), while in the case of twist-2 negative index 
sums were present starting from the second loop [13]. 

At the iV — > 00 limit, the minimal anomalous dimension 7 (corresponding to the 
ground state) must exhibit the universal (LBK-classical) In N behaviour which de- 
pends neither on the twist, nor on the nature of fields under consideration [25]. Com- 
puting analytically the large N expansion {1.7} yields (for /3 = Cs) 



7(iV) ~ ^4^\nN, 
Sph = r, -9 -Q29 + -5-C2 9 - I10C2 + j 5 + ■ ■ ■ , (2.iU) 



which matches the four-loop cusp anomalous dimension [22,26]. This is a non-trivial 
check, since the derivation of the expressions (|2.9|) was based on experimenting with 
finite values of the spin N. 
Expanding \\.2} gives 

7 = p + IVV + \ [PP^ + ^P^P] + I [PP^ + IV^PP + iP^V] + O {g^°) (2.11a) 
and, equivalently, 

7^ = 7 - 177 + 1 [7f + h'7] - -s [77' + i7'77 + b'''!'] + O {g'') , (2.11b) 
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where each dot marks derivative over N. 

Substituting 7 and V in terms of perturbative series in = 



n=l 



after a short calculation one obtains 

4 5i, 



p(i) 

p(2) 
p(3) 



-2 53 - 4C2 5i, 

55 + 2 C2 ^3 + 4 (53,2 + 54,1 - 2 53,1,1) 
+ 4 5i (2 53,1 - 54 + 4 C4) - 4 Sf (53 - Cs). 



(2.12a) 
(2.12b) 

(2.12c) 



The fourth loop kernel we split into two terms, the first built up of harmonic functions 
S only, and the second one containing the zeta-function factors. 



p(4) 



+ P, 



(4) 



C 



(2.12d) 



Here we present the corresponding expressions in terms of nested harmonic sums 
with unity indices moved to the head of the index vector, which form turns out to be 
more compact and is better suited for further transformation: 



,(4) 



.(4) 



— 8 [ — (53,3 + 5i,5 + 252,4) + 4(52,1,3 + 5'i,2,3 + 'S'1,1,4) — 85i, 1,1,3] 5i 

+ I 57 - 16 (5i,6 + 54,3) - 24 (52,5 + ^3,4) 

+ 48 (5i,i,5 + 5i,3,3 + 53,1,3) + 64 (52,2,3 + ^2,1,4 + 5i,2,4) (2.12e) 
- 128 (5i, 1,1,4 + 52,1,1,3 + ^1,2,1,3 + 'S'l, 1,2,3) + 256 5i, 1,1, 1,3; 

= 8C4 Sf - 4 [C2C3 + 8C5] ^i' - [4(C3 + 2;0)<S3 + 49C6] <Si 

+ (8<Si,i,3 - 45i,4 - 4^2,3 - Ss) C2 - SSs Ci- (2.12f) 



2.2 RR harmonic functions 



In order to show that the twist-3 evolution kernel (|2.12|) satisfies the Gribov-Lipatov 
reciprocity, we need to introduce the corresponding basis of harmonic functions. 
The functions 



-In^^x = r{2k + 1) ■ S2k+iix) 

X — 1 

satisfy the GL reciprocity il.5} and generate harmonic functions with an odd index. 



r(2A; + 1) 
1 



r{2k + 1) 



M 
M 



^In^ 

, a; — 1 

X 



S2k+i{N)] (2.13a) 
S2k+i{N) = S2k+i{N) - <S2fc+i(oo). (2.13b) 
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In ( |2.13bl ) and hereafter the hat marks the function in the Mellin moment space with 
its value at iV = oo subtracted: a{N) = a{N) — a(oo). 

The inverse Mellin transforms <Sni(a;) of multi-index sums SmiN) that enter eqs. 
(|2.12c|l - (|2.12f|l mix under the reciprocity operation (|1.5|) , so we need to construct linear 
combinations of harmonic fimctions <S that do have definite GL parity. 

2.2.1 Integral representation 

An index vector m = {mi , m2 , . . . , m^} corresponds to transcedentality r = Yli=-]. I^i I • 
We will refer to the transcedentality minus the length of the sum, w [m] = r — ^, as the 
weight of the harmonic function. 
Consider the recurrence relation 

hm{x) = -im]-'^ C ^^^i^ . *^(z), (2.14a) 

where the single index function coincides with the image of the standard harmonic 
sum, 

^aix) = [r(a)]-^ In"^-^ - = Sa{x). (2.14b) 

a; — 1 X 

At the base of the recursion, (|2.14b|) , we have 

^a{x) = l^-X^a{x-')y{-ir-'=(^-X^a{x-')y{-ir^''l 

An iteration (|2.14al ) increases the transcedentality of the function by b, and the length 
of the index vector by one, so that 

w[in] + b — 1 = w[b,m]. 

Observing that the integration measure in (|2.14a|) transforms under z — ^ 1/z as 

d^(z) = ^^ilJliZ^ #(z-^) = -z #(z), 

we conclude, by induction, that the functions ^ra{x) that the equation (|2.14|) generates, 
have definite GL parity determined by their weight: 

~^^{x) = (-l)-H .(^-x ^m{x-')j . (2.15) 

Let us represent the r.h.s. of (|2.14a|) as 

1 X f'^ dz . . z z - 1 ~ , . 2 X dz , f, z ~ , . ,^ , ^, 
I [bj X — 1 Jx z X z L [bj X — 1 Jx z X 
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Acting on a complementary harmonic su n the first integral produces Sa+b,n 

while the second generates (—2 times) the sum •Sa,b,n, see (|A.7|) . Iterating (|2.14|) , 

*c,d = Sc+d - ^S^d, (2.17a) 

*6,c,d = Sb+c+d - 2 { Sb+c,d + -^fc.c+d ) + ^ ^b^cd , (2.17b) 
^a,b,c,d = Sg+b+c+d — 2 (<Sa+6+c,d + Sa+b,c+d + <5a,6+c+£i) 

+ 4 { Sa-^-b^c,d + '^a,6+c,d + '^a,6,c+d ) " 8 Sa^b,c,d , etC. (2.17c) 

In the Mellin moment space, the formal construction of functions ^m{^) in terms of 
complementary harmonic sums is described in Appendix, see (|A.1H) , (|A.12|) . 

Equation (|2.15|) shows that the functions l>in with even weight vj[m] = t — I are 
Reciprocity Respecting (RR). 



2.3 Answer 

In terms of the functions ^ so introduced the evolution kernel (|2.12|) reads 

p(^) = 4<Si, (2.18a) 

P(2) = -4C2<Si-253, (2.18b) 

P(3) = ^^2_5^^2C2<S3 + 3,S5-2<i>i,i,3, (2.18c) 

p(4) ^ ^ifc,l + ^C,l)Si-C,2{^S^-2^,,i,^)-^e2S^ + ^SiM + B^, (2.18d) 
where 

M = 2 ii,i,i,3 - ($1,5 + $3,3) - i2p - Cs) S3, (2.19a) 

^4 = 16<i>l,l,l,l,3-4('l>3,l,3 + *l,3,3 + *l,l,5) -§57, (2.19b) 

We recall that the arguments of the functions on the r.h.s. of the equations is N/2, 
and P = in (|2.19a|l . Alternatively, in terms of the physical coupling (|2.10|) . the 
perturbative series for the kernel, V = J2n^i Sph ^ph'' takes the compact form 

r^i! = 45i, (2.20a) 

V^^ = -2 53, (2.20b) 

= 3Ss - 2*1,1,3 + C2 • (-253), (2.20c) 

= 4 5i • ^4 + + 2 C2 • (3 55 - 2 <l>i,i,3) . (2.20d) 

Since all harmonic functions involved have even weights w, the splitting kernel (|2.20|) 
is Reciprocity Respecting. Mark that given the argument N/2 of the harmonic func- 
tions in l|2.20|l , (|2.19|) , the reciprocity relation (|1.5|l applies to V{x^). 

This result can be compared with the evolution kernel that generates the twist-2 



See Appendix for the definition and properties of complementary harmonic sums 
*the GL parity of a harmonic function with k negative indices is (—1)'""'"'^ 
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universal anomalous dimension [15]: 



,(1) 
ph 



45i(iV)i 



(2.21a) 



,(2) 
ph 



(2.21b) 



V"^^ = 8 55(iV)-24*i,i,i,_2(^)-8C2<S3(iV) 

- 8 5i(iV)- [2*i,i,_2(iV) + $_2,-2(iV)-5_4(JV) + C2<S-2(i\r)]. (2.21c) 

2.4 Asymptotic expansion 



Large Mellin moments N correspond to parton light-cone momentum fractions a; — > 1. 
In this region only partons with small energy-momentum can be produced in the final 
state of a hard process, (1 — a;) ^ 1, and radiation of soft gluons dominates the answer. 
Therefore, the large N behaviour is important for understanding how the LBK physics 
manifests itself in the structure of anomalous dimensions and of the evolution kernel. 

At large positive values of the argument, z — > +oo, the harmonic function <Si(z) 
has the asymptotic expansion 



Si(z) = ^{z + 1) - ^(1) = (Iny + 7s) + \y-^ - ^y-' + ^y-' + ..., (2.22a) 



Such a structure of the series is inherited by all harmonic sums Sa{z), which can be 
obtained from (|2.22|) by simple differentiation. 



as well as by the nested harmonic functions "I'm (2). Namely, the asymptotic series of 
the functions of even (odd) weight contain only even (odd) inverse powers of y. This 
is a direct consequence of their inverse Mellin images having definite GL-parity [19]. 

2.4.1 Conformal structure of the large N expansion 

This feature of the asymptotic expansion has a clear symmetry origin [19]. 

On the light-cone the residual conformal invariance is that of the collinear sub- 
group of the conformal group 5L(2;R) C 50(2,4) [27]. In the j\A = 4 SYM theory, 
(super)conformal invariance is exact at all loops and quasipartonic operators can be 
classified according to the representations of 5L(2; M). Also, the operators belonging 
to different SL{2; M) multiplets do not mix under renormalization. The anomalous 



where 



y2 =z(z + 1). 



(2.22b) 
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dimension of each conformal multiplet depends only on the conformal spin charac- 
terizing the representation 



j{g) = i(iV + A{g)) = N + lL + \ jl{N), (2.23) 

where A{g) is the scaling dimension, A = N + L + Jl, which includes the quan- 
tum anomalous contribution "fhi^] g)- Hence, the spin dependence of the anomalous 
dimension at fixed twist L is encoded in the non-linear relation 

lL{N) = fL{N + \^L{N)], (2.24) 



which coincides with the relation (|1.2|) between the anomalous dimension and the 
Reciprocity Respecting evolution kernel V for the space-like evolution, a = —I. 

As we have mentioned in the Introduction, the fact that the asymptotic series for 
the RR evolution kernel runs in integer powers of (modulo logarithmic corrections) 
has been dubbed "parity preserving asymptotic series" and verified for a broad vari- 
ety of QFT models in [19] for the twist-2 case. Now we see that this property holds for 
twist-3 as well. 

In general, the conformal Casimir operator reads [28, 29] 

J'^ = {N + Lr}-1){N + Lr]), (2.25) 

where 77 = ^, 1, | for scalars, gaugino, and gauge boson fields, respectively. 

For twist-2, the argument of harmonic functions in eqs. \2.21) coincides with the 
Lorentz spin, z = N, and we have simply = N{N + 1) = y^. For the twist-3 in the 
scalar sector, (|2.20|) , we have z = N/2} substituting L = 3, 77 = | into (|2.25|) gives 

J^ = {N + i)(iV + I) = 4 (f + |)(f + I) = 4z2 + 4z + I = 4y2 + |. 

Thus, the expansion in y^ translates into that in the quadratic Casimir, J^. 

Asymptotic series for the evolution kernel, both for the twist-2 and twist-3 opera- 
tors, can be cast in the form 



P{N) = (Iny+jE) 



00 



n=l ^ 



bm \ y^ = N{N + l), L = 2; 



^-oV'"-' \y^ = f(f + l), L = 3. 

(2.26) 

where the L-dependent coefficients a, b are given in series of the coupling. 



2.4.2 Logarithmic structure of the large N expansion 

Apart from "parity preservation", the expansion (|2.26|) has another remarkable fea- 
ture. In a sharp contrast with the series for the anomalous dimension, where the number 
of logs increases with the power of the 1/N suppression, see (|1.7|) , the evolution kernel 
appears to be linear in In N. 
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Whether this property holds beyond the known orders {g^ for L = 2; for L = 3) 
remains unknown. Should this be the case, that would mean that the tower of sub- 
leading logarithmic singularities in the anomalous dimension are actually inherited, 
in all orders, from the first loop — the LBK-classical 7'^^ = V^^^ oc S\, and the RR 
evolution equation (|1.2|l generates them automatically. 

The leading logarithmic behaviour is universal and is given, in terms of the phys- 
ical coupling l|2.10|l , by the LBK-classical contribution to the evolution kernel, 

V{N) ~ 4g2j,5i(z) =4g2i^lni\r + (2.27) 

The non-trivial JV^ dependence accompanying the log N enhancement, 5V oc log N/N"^, 
appears only at the level of ai oc g^ and ai a g^ for twist-2 and 3, correspondingly. 

It follows from the structure of the relation (|2.11a|) , 

7W = E^ 5:4 lnN)]\ (2.28) 



, k\\2dN 

k=l 

that the maximally logarithmically enhanced contributions to the anomalous dimen- 
sion, log*^ N/N^, are driven by the leading, first loop term (|2.27|) . Representing the 
twist-3 anomalous dimension as 

M 00 . fc 

7(i\r;L = 3)=ao(ln-+7s)+S + ^-^ ^afc,^ln-i\r, ao=Agl^, (2.29) 

fc=l m=0 

and substituting (|2.27| > into ( I2.28D we immediately obtain 

1^2 „ i„3 „ 1 „4 . „ ("I)*' ^ „k+l CO on^ 

^1,1 — 2^0 ' ^^2,2 — - 8^0 ' ^^3,3 — 24*^0 i ^k,k — ^^"^ O-Q ■ (Z.^Uj 

In order to predict the next-to-maximal logarithmic terms, ak,k-i log*'^'^ N/N^, in the 
expansion (|2.29|) , it suffices to keep the constant in (|2.27|l and use the approximation 



V{N) ::,Ag'p^S,{z) + B = Ag^^- (In ^z{z + 1)+Jb] +B + (z- 

+ B, 



lnz + 7^ + ^ 



(2.31a) 



where 

B = S(5^ L = 3) = -2C3 gjh - (2C2C3 + Cb) g?h + O (g«h) . (2.31b) 
In particular, this gives (recall that for twist-3 we have z = N/2) 

ai,o = aoil + lB). (2.32) 

Relations l|2.30|) , (|2.32|) hold in all orders in the coupling constant. 
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3 Conclusions 



The notion of the reciprocity respecting evolution equation (RREE) emerged in an at- 
tempt to combine, in a single framework, anomalous dimensions of space-like parton 
distributions and time-like parton fragmentation functions, in order to simplify the 
structure of the higher order corrections. The basic observation is that the complex- 
ity of higher loop contributions is, to a large extent, inherited from lower orders. This 
is especially so for the major part of the QCD anomalous dimensions which is gov- 
erned by the "classical" gluon radiation, in the sense of the Low-Burnett-KroU (LBK) 
theorem. One may argue that there should exist a framework in which effects of the 
classical gluon fields would be fully generated, in all orders, from the first loop, in the 
spirit of the LBK wisdom. 

In the space of Mellin moments of parton distributions, the corresponding part 
of the anomalous dimension, in the nth order of the perturbative loop expansion, is 
described by harmonic functions of the maximal transcedentality , r = 2 n — 1. 

Since QCD shares the gauge boson sector with supersymmetric theories, the latter 
provide a perfect ground for testing the inheritance idea. The maximally supersym- 
metric JV = 4 SYM model is exceptional in this respect. In this theory quantum effects 
due to gluons, gauginos and scalars cancel in the beta-function (in all orders) as well 
as in the one-loop anomalous dimension. 

In this paper we have addressed the question of applicability of the notion of RREE 
beyond the leading twist. We considered the bosonic sl{2) sector of the jV = 4 SYM and 
found that the Gribov-Lipatov reciprocity relation holds in four loops for the evolu- 
tion kernel V describing the minimal (ground state) anomalous dimension of twist-3 
operators built of three scalar fields. The GL relation holds with substituted for 
X, in accord with the fact that the harmonic functions entering the twist-3 anomalous 
dimension have N / 2 as the argument. 

The twist-3 evolution kernel is given in eqs. (|2.20|) , (|2.19l ). Having derived this 
kernel from the space-like anomalous dimension, prediction for the time-like twist-3 
anomalous dimension follows immediately from the RREE. It is given by (|2.11al > with 
the r.h.s. changed to — r.h.s.[— P] (i.e., changing signs of the terms even in V). 

The structure of the twist-3 evolution kernel turns out to be even simpler than 
that for the leading twist-2, (12.211) . In particular, an extreme simplicity of the second 
loop kernel (|2.20b|l is likely to be due to the fact that the genuine O (a^) graphs that 
renormalise twist-3 operators vanish identically because of the colour structure (each 
pair in the adjoint representation): 
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In the large N limit, where classical LBK gluon radiation dominates, the asymp- 
totic series for evolution kernels are much less singular than the corresponding series 
for the anomalous dimensions. Namely, only the first power of 5i a In iV is present in 
the evolution kernel, both for twist-3 (|2.20|) , and twist-2 (|2.21|) . Whether this remark- 
able property is of general nature and holds in higher orders (n > 4 for twist-2, and 
n > 5 for twist-3) remains unknown. 

In twist-3, the subleading logarithmic enhancement <Si • Ai{N) oc In N/N^ appears 
for the first time at the level of a^, while in the case of twist-2 such term is present in 
the third loop, ~ In N/N^. 

The specific features of the minimal twist-3 anomalous dimension, and of the cor- 
responding evolution kernel, such as the absence of negative indices and appearance 
of N/2 as the argument of harmonic structures are awaiting physical explanation. 

Acknowledgements 

One of us (YLD) is indebted to Gregory Korchemsky and Benjamin Basso for an illu- 
minating discussion. 

A RR harmonic structures 
A.l Complementary harmonic sums 

Let m denote a string of integers, m = {mi, m2, . . . , m^}. We define multi-index (or 
"nested") complementary harmonic sums cS^ by the recursive relation 

Sa,m{N) = Sa,mW " Sa{N) ■ S^{cX^), (A.la) 

where we set 

Sa{N) = Sa{N), (A.lb) 

For any index vector m with the last index ^ 1, the value <Sin(oo) is finite, so that 
we can carry out subtraction to introduce 

Sn,{N) = SniiN) - <Sjj,(00). (A.2) 

For positive integer N we have 

oo 

Sa,miN) = - (sgna)"n-l'^l4,(n). (A.3) 
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From now on we restrict ourselves to positive indice^ so that for the index vector 
m = {mi , m2 , . . . , m^}, with mi > 2, 

JV oo oo 

S^{N) = (-1)^-1 ^n^^ J2 ^2^^' ■■■ E (A.4a) 

ni = l n2=ni + l n^=n^_i+l 



oo oo oo 



ni=N+l n2=ni+l n£=n^_i+l 

A.1.1 Large N behaviour 

As follows from their definition ( |A.4|) , (subtracted) complementary sums fall in the 
N ^ oo limit. The leading power of the large N behaviour is given by the weight of 
the sum, w (equal transcedentality minus length): 

i 

Srn{N) o^ N-'"^"'^ ■u;[m] =r-^ = E(m, - 1). (A.5) 

t=i 

Importantly, large N asymptotic expansion of complementary sums contains no log N 
enhanced terms, at any order of the 1 /N suppression. 

A.1.2 Mellin image 

Inverse MeUin images of complementary sums, 

Srnix) = [ Srn_{N) ] , (A.6) 

can be generated by the integral operation 

Sa,m{^) = [r(a)]-i^ [' — In^-i - • <^(z). (A.7) 

In particular, 

4a(a;) = [r(a)r(6)]-i— - / (A.8a) 

X — \ Jx Z — \ X Z 

s,,,M = [r(a)r(6)r(c)]-^-^ /' m-^ f /' -^m^-^ ^m-^ i (A.sb) 

x — lJx y—i- X Jy z — 1 y z 

etc. 



^Generalisation is straightforward; for complementary sums with a negative index see [15]. 
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A.1.3 Derivative over N 



■^SrniN) = M \ \-nX ■ S^{X)^^ = - ^ mfc • ■Smi,.,mfc,m,+l,m,+i,.,m, (iV), (A.9) 

fc=l 

Using the definition (|A.la|) , it is straightforward to derive the formula for the deriva- 
tive of the standard nested harmonic sums, 

d ^ - ^ - 



dN 

k=i k=i 



(A.10) 



A.2 Construction 



Given a complementary sum .Sm of length I, we construct a combination of sums 
y^'^^ having the same transcedentality and reduced length. The index vector m has I 
indices, separated by ^ — 1 commas; by erasing [I — r) of them, and summing up the 
non-separated indices, we construct 

Vm ~ "^m = Smi,m2,m3,...,mi ) (A. 11a) 

ym ^ ~ ^[mi+m2],m3,...,me ~l~ ^■mi,[m2+m3],...,me + • • • + '^mi,m2,...,[m^-i+m^] i(A-llb) 
Vrn ^ ~ ^[mi+m2+m3],m4,...,mi + '^[mi +7712], [7713+7^4],. ..,m£ 

+'^[77ii+77i2],m3, [7714+7715],. ..,m^ -|- . . . -|- <Smj^m2,...,[m£_2+"^£-l+m£]' (A.llc) 



i-1 



y(2) 
I'm 



S,^k . (A.lld) 



y^A^ = 5^ = -^r. (A.lle) 



Now, we consider definite combinations of complementary harmonic sums: 

*m = E(-2)'-SW = {-2Y-'S^ + ... + Ayg) - 2y(^) + Sr. (A.12) 

r=l 

Specificity of the combinations so organised lies in the structure of their inverse Mellin 
images. Namely, inverse Mellin images of (subtracted) harmonic functions (|A.12|) , 

^m{x) = M-^ l^UN)] , (A.13) 

have definite GL parity, 

^mix) = i-ir-(^-X imix-')^ . (A.14) 

This shows that the functions I'm with even weight io[m] = t — I are RR. 
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A.2.1 

r = 5: 



Some useful examples 



*1,1,3 = 4<Si,i,3 - 2(5i,4 + 52,3) + S5. 



(A.15) 



3,3 
1,5 



1,1,1,3 



-2 53,3 + Sq, 



(A.16a) 
(A.16b) 



-2 cSi,5 + Sq; 

-8 51,1,1,3 + 4(5i,i,4 + 51,2,3 + 52,1,3) - 2(5i,5 + 52,4 + 53,3) + 56. (A.16c) 



*i,3,3 = 45i^-2(5i^ + 54^) + 57i (A.17a) 

*3,i,3 = 453^-2(54^ + 53^) + 57i (A.lTb) 

*i,i,5 = 45i^-2(5i^ + ^) + 57i (A.17c) 
*i, 1,1, 1,3 = 16 5i, 1,1, 1,3 - 8( 5i, 1,1,4 + 5i, 1,2,3 + 5i,2,i,3 + 52,1,1,3 ) 

+ 4( 5i,i,i,4 + 5i,i,2,3 + 5i,2,i,3 + 52,1,1,3 ) (A.17d) 

- 2(5i,6 + 52,5 + 53,4 + 54,3) + 57. 
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